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xTliis  technical  summary  report  was  prepared  by  the  Research  Department  of  the  Allison 
^  Division  of  CJeneral  Motors  Corporation,  The  work  reported  was  accomplished  under  Contract 
Nonr-4i 04(00)  and  Amendment  No.  1  thereto. 

The  program  w'as  sponsored  by  the  Advanced  Research  and  Project  Agency  through  the  Power 
-  Branch  of  the  Office  of  Naval  Research  under  the  direction  of  ARPA  and  Mr.  J.  A.  Satkowski 
of  ONR. 

Allison  personnel  who  contributed  material  for  this  report  are: 

H.  E.  Wilhelm 
R.  O.  Whitaker 
R.  L.  Koch 


1 


TABLE  OF  CONTENTS 


Ser;tion  Title  Page 


I  Introduction . 1 

11  Summary  of  Work  Accomplished .  3 

Construction  of  Closed  Loop  Device  .  . .  3 

Test  Runs  . . .  . .  3 

Theoretical  Investigations .  4 

III  Experimental  Results .  5 

Runs  with  Pure  Helium . 5 

Experimental  Technique . 5 

Data  with  Pure  Helium . . 8 

Runs  with  Helium  Seeded  with  Cesium .  11 

Voltage-Current  Curves . 11 

Spectroscopic  Observations . 15 

IV  Theoretical  Investigation  . .  21 

Calculation  of  Intercomponent  Nonequilibrium  in  a  Magnet oplasma- 

dynamic  Converter . 21 

Converter  Model .  21 

Plasma  Equations .  21 

Initial  Flow  Data .  22 

Results  and  Discussion . 23 

Further  Investigation .  26 

On  the  Nonobservability  of  the  Full  Open  Faraday  Load  Voltage .  26 

V  Conclusions . 29 


LIST  OF  ILLUSTRATIONS 


Figure  Title  Page 


1  Segmented  electrodes . 5 

2  Magnet  current-field  strength  calibration  curve .  6 

3  Instrumentation  for  field  strength  recording . 7 

4  Voltage  output  recording  with  500 -ohm  load  . .  8 

5  Voltage  output  recording  with  1000-ohm  load . . .  8 

6  Voltage  output  recording  with  5000-ohm  load . . 9 

7  Voltage  output  recording  with  20,  OGO-ohm  load . 10 

8  Voltage  output  recording  without  load .  10 

9  Voltage-current  curves  for  cesium-seeded  helium .  12 

10  Voltage  output  curves  obtained  with  segmented  electrodes.  . .  14 

11  Spectrum  photograph  showing  reversal  of  the  cesium  resosnance  line .  17 

12  Densitometer  traces  showing  reversal  of  the  cesium  resonance  lines .  19 

13  MPD  converter  cross  section  model . 21 

14  Results  of  numerical  integration .  24 


LIST  OF  TABLES 


Table  Title  Page 


1  Elementary  particle  constants . 22 

11  Initial  flow  values  . . 22 

III  Electrical  flow  quantities  .  .  .  ,  . . 26 


V 


L  INTRODUCTION 


lliis  report  describes  the  progress  achieved  on  Contract  ONR  Nonr-4104(00)  from  1  August 
through  31  October,  During  this  period  sufficient  operation  experience  was  gained  and  design 
information  was  accumulated  to  justify  building  a  large  MPD  test  section.  This  new  test  section 
should  deliver  a  nominal  amount  of  power. 

The  Allison  MPD  system  is  now  sufficiently  refined  that  relia;  •  operation  can  be  ensured  for 
an  extended  time.  The  next  quarterly  report  will  describe  the  oesign  of  the  new’  test  section. 
Details  of  voltage-current  characteristics  and  saturation  curves  are  included  herein. 
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II.  SUMMARY  OJ  WORK  ACCOMPLISHED 


a  CONSTRUCTION  OF  CLOSED  LOOP  DEVICE 

^  Minor  design  modifications  in  the  closed  loop  system  were  made  during  the  reporting  period. 

I  The  helium-cesium  separation  system  was  improved  so  that  it  is  now  possible  to  recover  a 
large  percentage  of  the  cesium. 

’  The  seed  injection  system  has  also  been  improved  so  that  a  controllable  cesium  mass  flow  with- 
=  out  fluctuations  can  be  established  and  maintained  for  an  extended  time. 

The  design  and  construction  of  the  new  test  section  with  a  channel  length  of  50  cm  was  continued. 

~  TEST  RUNS 

s  Several  helium-cesium  runs  were  made.  Voltage-current  curv^es  for  different  magnetic  fields 
were  plotted  directly  with  x-y  recorders.  A  family  of  curves  with  the  magnetic  field  as  a 
parameter  was  thus  obtained.  The  voltage-current  curves  are  straight  lines.  'Die  magnetic 
I  field  causes  the  voltage-current  curves  to  rotate  around  one  point  at  the  current  axis,  as  is  ex¬ 
pected  fr^ni  theory.  Direct  plotting  of  the  voltage-current  curves  ensures  that  the  curves  are 
^  straight  lines  and  have  no  occasional  dips  which  could  be  overlooked  in  the  interpolation  when 
composing  curves  from  different  runs. 

g  For  the  same  reason,  efforts  were  made  to  record  the  voltage  output  curves  versus  magnetic 
I  field  directly  from  the  generator  to  x-y  plotters.  Ln  this  case  the  load  re  istor  was  the  param- 
I  eter.  Since  it  is  not  possible  to  measure  the  magnetic  field  directly  during  the  run,  the  cur- 
I  rent  through  the  magnet  coils  was  monitored  and  this  value  was  converted  with  a  calibration 
I  curve  into  the  corresponding  value  for  magnetic  field.  For  a  direct  plot,  this  transformation 
I  must  be  done  continuously.  'Dierefore,  the  signal  from  the  magnet  current  shunt  was  sent  via 
J  telephone  lines  to  the  computing  center  where  an  analog  computer  converted  it  to  a  signal 
I  corresponding  to  the  magnetic  field.  This  signal  was  then  fed  back  into  the  x-y  plotter. 

Different  spectroscopic  measurements  were  also  made  and  a  few  impurities  were  detected. 

The  essential  result  of  the  spectroscopic  observations,  however,  is  that  the  electron  tempera- 
jture  was  equal  to  the  gas  tempera  ure  within  the  accuracy  of  measuement. 


TH  EOR  R  TIC  A  L  I NV  ESTIG  A  TIONS 

A  theory  of  the  intercomponent  thermal  nonequilibrium  was  developed  at  small  Hall  coefficients 
including  nonuniformitiec  of  the  plasma  and  the  thermal  diffusion  effect. 

First  numerical  results  for  the  buildup  of  the  electron  temperature  noriequilibrium  were  obtained 
using  the  usual  approximation  which  disregards  the  thermal  diffusion  effect.  The  numerical 
value  of  the  buildup  length  is  in  good  agreement  with  the  approximate  analytical  formula  ob' 
tained  earlier. 


RUNS  WITU  JRE  HELIUM 


A  previous  discussed  MPD  runs  made  with  a  workinjf  fluid  of  pure  helium  with 

Observations  made  during  runs  of  that  type  are  reported  herein. 


Experimental  Technique 


In  the  original  run,  one  of  tiie  electrodes  (the  cathodf  )  was  a  single  plate,  while  the  other  elec¬ 
trode  (the  anode)  consisted  of  pins.t  In  the  present  setup,  segmented  electrodes  were  used  and 
only  the  high-voltage  electrodes  were  pins.  Figure  1  shows  the  electrode.s. 


Figure  1.  Segmented  electrodes 


Measurement  of  Fluid  Proprrtic  s  for  Magnetoplasrnadynamic  Power  Generators  —  Fourth 


Quarterly  Technical  Summary  Report  (1  February  —  3(}  April  1(164),  Indianapolis.  Indiana 


All'son  Division,  GMC,  FDR  3861,  16  May  1964 
!-.id.  .  Figure  1 . 
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Segmented  electrodes  were  used  to  prevent  the  cathode  spot  from  moving  downstream  as  had 
happened  in  previous  experiments.  *  With  an  arrested  cathode  spot  as  is  achieved  with,  the  pin 
electrodes,  it  is  much  easier  to  judge  how  far  downstream  the  induced  conductivity  can  be 
extended.  It  should  also  be  noted  tiiat  'with  the  high-voltage  pin  electrode,  the  electron  emitting 
rfuciace  is  several  orders  of  magnitude  smaller  than  previously  and  the  area  is  more  precisely 
defined. 

During  system  operation,  it  is  vt  ry  desirable  to  monitor  the  magnetic  field  strength  in  the  test 
section.  However,  no  transducers  are  available  .hich  can  be  positioned  in  the  test  section  and 
can  endure  the  environ'  The  most  practical  alternative,  then,  is  to  measure  magnet  current 

The  calibration  curve  wn  m  Figure  2.  may  then  be  used  to  obtain  field  strength.  To  avoid 
the  necessity  of  hand  reducing  data,  the  instrumentation  setup  of  Figure  3  was  employed. 


Ciirve  u  for  risini:  current  in  magnet 


Magnet  current  “  amp 


Figure  2.  Magnet  current — field  strength  calibration  curve. 


^Ibid.  ,  Figure  5, 
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Cesium 


Figure  3.  Instrumentation  for  field  stren,^th  recording. 


Magnet  current  passing  through  the  50-mv  shunt  developed  a  voltage  which  was  amplified  by 
the  differential  d-c  amplifier.  ^  This  amplifier  is  a  chopper-stabilized  unit  of  very  low  drift 
and  very  constant  gain. 

A  specially  install^  •  cable  delivered  the  output  of  the  amplifier  to  the  Allison  analog  computer 
facility — indicated  by  the  portion  of  Figure  3  within  the  dotted  perimeter. 

The  gain  of  the  differential  d-c  amplifier  was  adjusted  so  that  300  amp  flowing  through  the 
50-mv  shunt  caused  the  output  of  amplifier  29  to  change  from  0  to  100  v.  Amplifier  TO  pro¬ 
vided  further  gam.  Potentiometer  42  (associated  with  amplifier  30)  compensated  for  any 
slight  zero  drift  in  the  differential  d-e  amplifie’'.  The  +i00-v  input  and  the  gain  of  2  applied 
to  the  main  signal  provided  proper  scaling  for  the  function  generator  56.  Function  generator 
56  carried  tiie  magnetic  field  strengtli  vs  current  curve  of  Figure  2.  The  output  of  the  function 
generator  was  thus  a  voltage  proportional  to  magnetic  field. 

The  output  of  function  generator  56  was  .  'turned  via  cable  to  the  MPD  laboratory  and  delivered 
to  the  x-axis  of  an  x-y  recorder.  This  permitted  variables  to  be  plotted  directly  as  functions 
of  magnetic  field. 


^Kin-Tel  Model  li4C  d-c  amplifier  and  Model  ^OlA  rack  module. 
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It  is  noted  that  the  magnet  used  exhibited  substantial  hysteresis.  The  curve  set  up  on  the  func¬ 
tion  generator  was  that  for  rising  magn^  ■  sation.  Proper  operation  of  the  instrumentation  re¬ 
quired  that  before  each  test  the  magnet  excitation  be  reduced  to  zero.  Magnet  current  was  then 
raised  to  approximately  the  test  value  and  held  for  the  duration  of  the  test. 


Data  with  Pure  Helium 

Figures  4,  5,  6,  7,  and  8  represent  the  raw'  data  obtained  from  low- temperature  runs  with 
helium  without  seeding.  The  output  voltage  measured  across  a  load  re&istor  was  plotted  di¬ 
rectly  as  a  function  of  magnetic  field.  The  load  resistance  appeared  as  a  parameter.  These 
curves  are  presented  here  because  a  saturation  effect  was  observed.  If  the  suspected  reason 
for  this  effect  is  substantiated  by  further  tests,  it  will  be  described  in  tiie  final  report.  At 
that  time,  runs  will  also  have  been  made  using  the  large  test  section  being  constructed  at  the 
present  time. 


Figtire  4.  Voltage  output  recording  with  500-ohm  load. 
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Figure  5.  Voltage  output  recording  with  1000-ohm  load. 
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Figure  6.  Voltage  output  recording  with  SOOO-ohm  load. 
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=  .RUNS  WITH  HELIUM  SEEDED  WITH  CESIUM 

iThe  analog  instrumentation  described  previously  was  used  to  obtain  Ixjth  voltage-  ci  ent  and 
“Voltage  output  curves. 

Voltage- Current  Curves 

Figure  9  shows  a  family  of  voltage -current  curves  taken  with  different  magnetic  field  strengths. 
The  curves  are  nearly  straight.  Except  for  very  low  magnetic  field  strengths,  the  magnetic 
field  seems  to  rotate  the  voltage- current  curves  around  one  point  on  the  current  axis.  This 
can  be  explained  as  follows. 

The  voltage  on  a  load  resistor  can  be  expressed  as; 

''l- 

where 


is  the  load  resistor 
Rj  is  the  internal  resistance 
V,  d,  £uid  B  are  constants 


The  current  is; 


I  =  [l/(R,  t  Rj^)] 


vBd 


(2) 


In  the  case  of  the  voltage -current  curves,  v.  d.  and  Rj  are  constant  and  B  is  a  parameter. 
Rj  may  be  a  function  of  B. 


Equation  (2)  can  also  be  written: 


Rj^  =  (vBd/1)  -  Rj 
and  therefore. 


(3) 


V,  ^  IR.  =  vBd  -  If: 

l-v  .w  i 


«) 


li 


iliif  .  '  i||i 


current  curves  for  cesium-seeded  helium 


Now,  assuming  r~  (as  was  explained  previously ' then  far  large  values  of  wr 

and  the  followang  expression  is  obtained; 


=  B(vd  -  C I  Rj,^) 


(5) 


whore 


Rjg  is  the  internal  resistance  without  magnetic  field 
C  is  a  constant 

Equation  (5)  show's  that  for  large  values  of  «t  ,  the  open  circuit  voltage  is  approximately  pro¬ 
portional  to  the  magnetic  field  w'hile  the  short  circuit  current  =  0)  does  not  depend  on  B. 
This  explains  the  observed  rotation  of  the  I-V  curves  around  one  point  at  the  short  circuit  cur¬ 
rent  axis. 


For  small  values  of  magnetic  field.  Equation  (5)  does  not  hold.  If  6>t  is  in  the  order  of  unitj% 
the  short  circuit  current  also  depends  slightly  on  the  magnetic  field,  which  can  be  easily  seen 

I  +  w  in  Equation  (4): 


I  =  vBd/Cg  Rjo  ^[l  ^  (e^B^  r^/m^) 


(6) 


w'here  Cg  =  constant 


Voltage  Output  Curves 


Figure  10  show's  a  family  of  voltage  output  curves  for  segmented  electrodes,  directly  plotted 
with  an  x-y  plotter.  The  character  of  the  curves  agrees  with  the  curves  previously  reported.^ 
Th  ’  curve.s  were  taken  w'ith  different  load  resistors  to  contrast  with  the  previous  data  where 
only  one  load  resistor  was  used.  To  compare  the  voltage-current  and  the  voltage  output  curves, 
lines  of  constant  resistance  were  drawn  on  the  voltage -current  curves.  Transferring  these 
points  to  the  voltage  output  curves,  on  a  line  of  constant  field  strength,  showed  that  the  points 
are  very  close  to  the  actual  recording. 


^'Measurement  of  Fluid  Properties  for  Magnetoplasmadynamic  Power  Generators  —  Fifth 
Quarterly  Technical  Summary  Report  (1  May  —  31  July  1964),  Indianapolis,  Indiana,  Allison 
Divi.sion,  GMC,  EDR  3965,  17  August  1964,  p  31. 
tlbid. ,  p  30, 
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Test  No. 


age  output  curves  obtained  with  segmented  electrodes 


^4  possible  explanation  for  the  saturation  effect  was  given  previously.  ^  After  runs  have  been 
:ffladc  using  the  large  test  section  now  under  construction,  more  information  will  be  available 
^to  support  or  uiSprove  this  explanation. 

SPECTROSCOPIC  Ob..ERVATIONS 

The  primary  problem  associated  with  spectroscopic  observations  is  the  very  weak  intensity  of 
^the  radiation  of  the  helium -cesium  mixture  at  approximately  ISOO'^K.  The  intensities  of  the 
lines  have  been  calculated  and  tabulated.^  As  has  been  demonstrated  previously,^  the  radiation 
of  all  lines  is  extremely  weak  at  the  temperature  range  under  consideration.  However,  it  w'as 
possible  to  detect  a  fairly  large  number  of  iines,  within  reasonable  exposure  times.  The  ex- 
^posure  time  for  tlie  resonance  lines  was  reduced  to  5  sec.  The  exposure  time  for  medium  lines 
is  approximately  1  min;  weak  lines  require  exposure  times  of  30  min  or  more. 

The  operation  of  the  closed  loop  system  is  now  reliable  enough  that  constant  gas  temperature 
-  and  constant  cesium  mass  flow  can  be  guaranteed  much  longer  than  30  min.  Those  cesium  and 
=  impurity  lines  which  have  been  detected  thus  far  are  as  follows- 


Cesium 

Rubidium 

Sodium 

8761A 

794  7 A 

8195A 

8521A 

780oA 

8 184  A 

8078A 

589oA 

7944A 

6973A 

4593A 

O 

4555A 

To  measure  the  absolute  value  of  the  electron  temperature,  both  the  line  reversal  metiiod  and 
the  line  ratio  method  have  been  used.  To  measure  relative  changes  in  electron  temperature, 
the  comparison  of  the  intensity  of  the  same  line  at  different  temperatures  proved  quite  sensi¬ 
tive.  This,  of  course,  can  be  seen  from  the  steep  slopes  of  the  curves  showm  previously. 


*lbid. ,  p  3 1 . 

fMeasur^ment  of  Fluid  Properties  for  Magnetoplasmadynamic  Power  Generators — First 
Quarterly  Technical  Summary^  Report  (1  May — 31  July  1963),  Indianapolis,  Indiana.  Allison 
Division,  GMC,  FDR  3511,  20  August  1963,  pp  45-58. 
t Ibid.  ,  Figure  19. 

"*Ibid. 
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The  line  reversal  method  is  well  known  and  is  often  .'Sed  especially  when  applied  to  the  reso¬ 
nance  lines  of  the  alkali  metals.  The  principal  shortcoming  of  the  line  reversal  method  lies 
in  the  assumption  that  the  temperature  distribution  is  uniform  within  the  plasma  volume.  If 
this  is  not  true^  then  the  errors  introduced  are  difficult  to  correct. 

Also,  extreme  care  must  be  used  when  aligning  the  optical  light  path.  •  Losses  due  to  absorption 
in  the  windows  and  in  Colder  gas  regions  in  the  light  path  may  introduce  errors.  This  method 
is  quite  sensitive  to  changes  in  temperature.  It  was  possible  to  differentiate  between  two  tem¬ 
peratures  less  than  50°K  apart  at  the  1200°K  level. 

The  absolute  value  of  the  temperature  is  relatively  inaccurate,  of  course,  because  of  the  error 
introduced  by  losses  due  to  absorption  and  the  error  in  measuring  the  comparison  blackbodv 
temperature.  By  using  very  conservative  estimates,  however,  the  error  in  the  temperature 
measured  with  this  method  cannot  possibly  be  larger  than  ±  T5°K. 

Because  of  the  question  of  temperature  nonhomogeneity,  line  reversal  was  recorded  photo¬ 
graphically.  Since  temperature  nordiomogeneity  can  now  largely  be  disconnected,  however, 
photoelectric  recording  is  used,  which  makes  routine  measurements  much  easier  and  faster. 

Figure  1 1  is  a  photograph  of  the  reversal  of  the  cesium  resonance  line  approaching  the  reversal 
temperature  at  IISOTC.  Figure  12  shows  the  densitometer  traces  of  these  photographs. 

As  a  second  method,  the  line  ratio  method  was  used.  If  a  temperature  profile  exists,  and  if 
the  plasma  has  cylindrical  geom.etry,  this  method  can  be  corrected  by  the  Abel  inversion. 

Here  the  line  pair  X  -  '1555A  and  X  -  807SA  v/as  used.  The  results  agreed  v/ith  the  results  of 
the  line  reversal  method  w'ithin  10%.  Of  course,  it  is  most  desirable  to  use  more  than  one 
line  pair.  To  do  this,  the  exposure  times  must  be  increased  further.  Runs  of  this  type  will 
be  made  in  the  near  future. 


Summarizing  the  spectroscopic  temperature  measurements  to  date,  the  electron  temperature 
was  always  equal  to  the  gas  temperature,  within  the  accuracy  of  measurement. 


16 


Fipure  11.  Spectrum  photograph  showing 
reversal  of  the  cesium  resonance  line. 


Fiij2.ure  12.  Densitometer  traces  showing  reversal  of  the  cesium  resosiance  lines. 
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IV.  THEORETICAL  INVESTIGATION 


CALCULATION  OF  L^TER COMPONENT  NONEQUILIBRIUM  IN  A  MAGNETOPLASMADYNAMIC 
CONVERTER 


Converter  Model 


The  converter  channel  is  assumed  to  be  a  rectangular  duct  of  constant  cross  section  as  shown 
in  Figure  13.  The  two  horizontal  duct  walls  are  solid  conductors  and  serve  as  electrodes.  The 
electrodes  are  connected  with  a  load,  R.  The  magnetic  field  of  constant  induction,  P^,  is  ap* 
plied  perpendicular  to  the  channel  axis  and  coplanar  to  the  electrode  planes.  The  working  fluid 
consists  of  four  components^ — helium  as  carrier  gas  (o),  cesium  as  seeding  gas  (a),  and  the 
ionization  products  of  the  latter,  single  cesium  ions  (i)  and  electrons  (e). 

Plasma  Equations 

For  the  theoretical  description  of  the  four -component  flov.,  the  multifluid  equations  for  reacting 
plasmas  based  on  the  13 -moment  approximation  are  used.  These  field  equations  have  been 
given  in  detail^  and  are  not  reviewed  herein.  To  reduce  the  numerical  effort,  a  few  simplifi¬ 
cations  are  made.  Viscous  effect.?  and  boundary  layers  are  disregarded.  The  Hall  effect  is 
disregarded  and  the  flow  is  treated  as  one -dimensional — i.e.  ,  all  fields  depend  only  on  the 
coordinate  in  the  direction  of  the  channel  axis.  (Assume  that  the  Hall  effect  is  compensated 


Figure  13.  MPD  converter  cross  section  model. 
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by  some  method.  )  FurtherTnorej  energy  losses  due  to  heat  conduction  are  not  taken  into  account 
Full  corisideration  is  given  to  <e-a)  ionisation  and  (e-i-e)  recombination  reactions,  effects  asso¬ 
ciated  with  the  n^^  ^uniformity  of  the  components,  momentum  transfer  by  intercomponent  friction, 
and  plasma  heating  by  intercomponent  friction. 

Initial  Flow  Data 


Two  converter  flows,  I  and  II.  vvith  different  exterior  loads,  R,  arc  considered.  The  exterior 
magnetic  field,  Bg,  is  the  same  In  both  cases.  Furthermore,  at  the  channel  entrance  (x  =  0), 
the  plasma  pressure.  P,  temperature,  T,  Mach  number,  M,  and  particle  composition  are  the 
.same.  The  ionization  at  the  channel  entrance  is  calculated  in  accordance  with  tlie  Saha  equation. 

The  strong  temperature -dependent  Gvosdover  cross  sections  (Qg^  =  are  calculated  in  the 

program.  The  cross  sections  for  neutral -neutral  iQ^o  ”  ^oa^  neutral-charged  (QgQ  =  Qgg 
-fid  Qga  =  Qgg,  where  s  -  e,  i)  particle  interactions  trigid  sphere  model)  are  shown  with  ele¬ 
mentary  particle  constants  in  Table  I, 


The  characteristic  quantities  of  the  fIov%f  at  the  channel  entrance  are  listed  in  Table  II. 

Table  I. 

Elementary  particle  constants. 


Q  =  5.  5  X  10'^^  cm^ 

=  4  X  iU  cm 

-  -  1  fi  2 
_  =  M  X  iU  cm 

Qjg  ^  2  X  icr^"^  cm^ 

(  ^  =  3 .  893  ev 


m^  =  6.  540  X  10  gm 
=  2.  172  X  10'^^  gm 
=  2.  172  X  10'^^  grn 
m^  -  9.  108  X  lO'^®  gm 
e  =  4.  803  X  10' esu 


Table  II. 

Initial  flow  values, 


6  2 

Plasma  pressure.  P  -  1,  5  x  10  dynes /cm^” 

Flow  velocity.,  v  =  1,  847  X  lU^  cm/sec 

Magnetic  induction,  B_  1=  18,  000  gauss 

^  1 H  - 1 

Helium  density,  n.  =  .5.38  X  10  cm 


Electron  density,  rr 


4.  24  X  10^'^  cm  ^ 


Flow^  I— R  2.  0  ohm 


Plasma  temperature,  T_  =  200n®K 

Mach  number,  M  0  8 

-2 

Seeding  ratio,  a  -  in 

Cesium  density,  r  -  5,  38  X  cm 

“  1  ■*  -  3 

Ion  density,  n  ^  4.  24  X  10  cm 

"  i 

Flow  11  —  R  -  6  6  ohm 


*S,  C,  Brown,  Basic  Data  of  Plasma  Physics,  J.  Wiley  and  Sons.  19a9. 
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Results  and  Discussion 

I  With  respect  to  converter  constants,  the  flows  I  and  II  differ  only  in  the  magnitude  of  the  load 
(Rj  <  Rjj).  This  has,  however,  considerable  consequences  for  the  flow'  properties  above  the 
^  channel  entrance,  w'here  x  >  0.  The  load  voltage  is  given  by  U  =  [R/(Rj+  R)]  vBd,  w'here  R, 

^  is  the  irRf^rior  resistance  of  the  duct.  Consequently,  the  field  impressed  by  the  load  into  the 

plasma  (E  =  -U/d)  is  smaller  for  flow  1  than  for  flow  II  w'ith  respect  to  the  modulus  (i  Ejl<iEjj|  ). 
~  Since  the  effective  field  acting  on  the  electron  compenent  is 

=E+^XB^d  (7) 

It  is  to  be  expected  that  flow  1  would  exhibit  stronger  nonequilibrium  effects  than  flow  II. 

Figure  14  shows  the  main  results  of  the  numerical  integration — the  dependence  of  the  electron 
temperature  (T^)  and  the  electron  density  (n^)  on  the  flow  coordinate  (x).  For  comparison  pur¬ 
poses,  the  gas  temperature  field,  T^,  is  also  gWen.  The  remaining  fields  of  the  heav'y  pa:  tide 
components,  which  did  not  change  appreciably,  are  not  shown.  In  :  ccordance  with  the  qualita- 
'  tive  considerations,  the  exact  numerical  calculations  indicate  a  more  than  twice  as  strong  bher- 
mal  nonequilibrium.  |  T^  -  T^l  ^  0.  for  flow  1  than  for  flow  11.  From  Figure  14,  the  charac¬ 
teristic  length  after  which  the  electron  temperature  is  built  up  is  read  as; 

Flow  I:  =  1.  1  X  lO'^  cm 

Flow  II  E'p  =  1-4  X  10  ^  cm 

According  to  a  previous  analytical  investigation,  1  is  approximately  given  by: 


=  v/(l6/3)  ^  2  k  Tp/ w  mg  {m^/mj,)(nj. Qj.) 


Because  of  the  predominating  (e-o)  interaction,  it  follows,  using  the  constants  from  Tables  1 
and  11,  that: 

irj,  =  v/(16/3)'^(2  k  T^/ir  m^Km^/ m^Rn^  (9) 

-  (2.  2)(10’2)  cm 


^assuming  <<  n,^  and  where  (r  =  o,  a,  i) 
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Figure  14.  Results  of  numerical  Integration 


As  th(>  electron  temperature  changes  conaiderablv  as  x  varies  between  zero  and  I it  is 

^  "  A  A  1 

plausible  to  insert  a  mean  value  defined  by  ^  (1/iHT.^  +  Tp),  where  T_  is  the  electron  tem¬ 
perature  at  X  into  Equation  (9),  With  this  correction,  the  agreement  with  the  numf-ricai 

calculations  becomes  fairly  good. 

From  Figure  14,  the  full  nonequilibrium  electron  temperature  ts  read  as 


Flow  1;  =  2846“K 

Flow  II:  Tp  =  2325'’K 


Compare  these  values  again  with  the  simple  analytical  theory  '  according  to  which 


Te/Tp  =  I  +  lEfV 


The  critical  plasma  field  for  frictional  electron  heatmg  is  given  by; 


w'hcre  the  electron  collision  frequency  (  v^)  is: 

•'e  =  <8/3)V2kTp/ir  (12) 

Because  of  the  predominating  (e-o>  interaction,  becomes,  with  the  constants  from  Tabic  II; 

Ef  S  (8/e)V(2/3ir)(mg/mo)<Tp/Tp)  Qeo^okT^j  (13) 

=  22.37-^Tp/T^  v7cm 

According  to  Equation  (10),  the  field  S-^  required  for  the  production  of  an  electron  temperature 
Te  =  2846’K  and  T^^  =  2325'’K,  respectively,  is  -  14.  53-^Tp/Tp  and  9.  (;5^Tp/T„  v/cm. 

respectively.  These  values  again  agree  satisfactorily  with  the  results  obtained  numerically. 


-Measurement  of  Fluid  Properties  for  Magnetoplasmadynamic  Power  Generators,  Fourth 
Quarterly  Technical  Summary  Report  (1  February  —  30  Aprii  1964),  Indianapolis.  Indiana. 
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Figure  14  does  not  indicate  a  noticeable  change  in  electron  density,  n^.  It  can  be  shown  that 
the  characteristic  length,  required  to  build  up  the  electron  density  corresponding  to  the 
elevated  electron  temperature,  is  large  compared  to  the  channel  length.  The  low-  electron 
density  explains  the  low  values  of  tlie  generated  power  for  the  two  flow  examples  (see  Table  111), 

Table  III. 

Electrical  flow  quantities, 

Loadj  R  (ohm)  Current  (amp)  Power  (w) 


Flow!  2.0  4.90  48.98 

Flow  II  6.6  3.03  60.53 

Further  Investigation 


These  preliminary  calculations  will  be  followed  by  computation  for  larger  converter  models  for 
which  the  build-up  of  the  electron  density  can  be  observed  and  a  noticeable  aeceleration  of  the 
flow.  Furthermore,  it  is  intended  tc  include  additional  physical  effects  which  may  be  of  im¬ 
portance  for  the  nonequilibrium  aspects  of  converter  flows. 


ON  THE  NONOBSERVABILITY  OF  THE  FULL  OPEN  FARADAY  CIRCUIT  VOLTAGE 

The  open  load  voltage  of  an  ideal  Faraday  converter  is  =  J" vB^dx^Iv  =  flow  velocity.  B^  = 

0 

magnetic  induction,  d  =  electrode  distance).  In  actual  experiment,  in  general,  an  open  load 

a 

voltage  Uj^  <  J'  vB^dxj^  is  observed.  In  the  following,  it  is  shown  that  this  can  be  explained 
0 

by  Hall  currents  closed  by  leakage  resistances.  Other  influential  effects  may  still  exist,  how¬ 
ever. 


For  Hall  coefficients  =  (tr /n^c)  Bq  ~  1,  the  Hall  effect  is  no  longer  negligible  in  plasma. 
In  this  case,  Ohm's  law  assumes  the  form; 


j  ♦  «T  j  X  (Bq/Bo)  -  <r  (E  +  V  X  Bq) 


(14) 
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It  follows  for  the  current  density  parallel  to  the  flovi  C  n  )  and  for  the  current  density  perpen¬ 
dicular  to  the  flow  and  magnetic  field  (i-): 


32  ^  io-/(l  +  w^T*-)]  [ej^  +7XB^-  «t  E,,  X<Bq/B^)] 


(15) 

(16) 


To  measure  the  open  load  voltage  of  the  Faraday  converter,  the  Faraday  load,  ,  is  removed 
and  replaced  b>'  a  high  impedance  voltmeter.  According  to  Equation  (16),  the  instrument  should 
indicate  the  voltage  (jj^  ^  0): 


d 


0 


whe  re 

-e2  =vXB^-  mr  E,,  X  (Bg/B„)  (18) 

The  field  E||  is  impressed  on  the  plasma  by  a  leakage  resistance. 

After  elimination  of  E,,  by  means  of  Equation  (15),  Equation  (18)  becomes; 

-e2  =  vX  -  (  4.r  X  (B^/B^)  (19) 

Equations  (Id)  and  (19)  show  that  in  the  presence  of  the  Hall  effect,  the  ideal  Faraday  field, 
-E^  -  "v  X  B,  is  diminished  by  a  field  -Ej^  -  -  ur  E,,  X  (Bg/B^)  -  jv)  X 

which  15  induced  by  the  flow  of  the  Hall  current  across  the  magnetic  field.  Substitution 
Equation  (19)  into  Equation  (17)  finally  gives  the  open  Faraday  load  voltage  in  the  form  (b  = 
duct  width); 

-(  «T  /erb)  I,,  (2U) 


where 

Iceal  open-load  voltage: 


Hall  current;  l|i 


(21) 

(22) 
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Because  «r/<r  s  B^/n^e,  the  observed  open  load  voltage,  U^.  should  decrease  as  the  electron 
density,  n^,  decreases  (at  low  temperatures)  and  increase  as  the  magnetic  field,  increases. 
To  experimentally  prove  Equation  (19),  the  Hall  leakage  current,  1,^  ,  must  be  measured.  A 
comparison  of  this  preliminary  theory  with  experiment  will  be  presented  in  the  final  report. 


V.  CONCLUSIONS 


SYSTEM  DESIGN 


The  major  problems  in  the  design  of  a  magnetoplasmadynamic  system  have  been  those  associ- 
ted  with  the  gas  heater,  cesium  seeding,  helium  purification,  and  helium-cesium  sep^i.  ation. 
These  problems  have  now  been  solved  to  the  degree  that  reliable  operation  of  tlie  MPD  system 
has  been  achieved. 

GENERATOR  DESIGN 

Operation  of  the  present  small  MPD  system  yielded  sufficient  information  to  permit  design  of 
the  large  test  section  which  is  now  under  construction.  Ceramic  and  tantalum  parts  were 
never  subject  to  severe  attack  by  cesium.  The  quartz  window  in  contact  with  tlie  hot  cesium - 
helium  mixture  was  found  to  devitrify  immediately.  No  attack  occurred,  however,  when  the 
windows  were  sufficiently  cooled.  It  was  necessary,  therefore,  to  remove  the  windows  from 
direct  contact  wRh  the  hot  working  fluid. 

ELECTRON  HEATING  EFFECT 

As  previously  m.entioned,  no  increase  in  electron  density  is  expected  in  the  4-cn.  channel 
(sec  Figure  14),  Therefore,  the  usefulness  of  the  electron  heating  effect  can  not  be  deter¬ 
mined  before  the  50 -cm  channel  is  in  operation. 

As  preliminary  indications,  however,  the  following  observations  can  be  reported.  It  was  not 
possible  to  observe  an  increased  electron  temperature  with  load  resistance  down  to  100  ohm. 
Since  the  measurements  reported  in  previous  quarterly  summary  reports  w'ere  made,  the 
numerical  results  represented  in  Figure  14  became  available.  From  these  results  it  can  be 
seen  that  an  increased  electron  temperature  can  hardly  be  expected  when  using  a  load  even 
as  low  as  100  ohm.  Attempts  are  being  made  now'  to  detect  this  effect  when  running  with 
short-circuit  loads. 
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